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ABSTRACT 

Using the Very Large Array and the Pie Town antenna, we have measured the position of the radio pulsar 
B 195 1+32 relative to nearby background radio sources at four epochs between 1989 and 2000. These data show 
a clear motion for the pulsar of 25 ± 4 milliarcsec yr" 1 at a position angle 252° ± 7° (north through east), corre- 
sponding to a transverse velocity 240 ± 40 km s" 1 for a distance to the source of 2 kpc. The measured direction of 
motion confirms that the pulsar is moving away from the center of its associated supernova remnant, the first time 
that such a result has been demonstrated. Independent of assumptions made about the pulsar birth-place, we show 
that the measured proper motion implies an age for the pulsar of 64 ± 1 8 kyr, somewhat less than its characteristic 
age of 107 kyr. This discrepancy can be explained if the initial spin period of the pulsar was Pq = 27 ± 6 ms. 

Subject headings: ISM: individual (CTB 80) — pulsars: individual (B1951+32) — radio continuum: ISM — 
stars: neutron — supernova remnants 



1. INTRODUCTION 

Associations between pulsars and supernova remnants 
(SNRs) allow measurements which would not be possible on 
either population of object alone. Since the center of an SNR 
marks the presumed pulsar birth-site, the pulsar's characteristic 
age, combined with the offset of the pulsar from the center of 
the SNR, lets us estimat e the pulsar's transverse velocity (F rail, 
Goss, & Whiteoak 1994). A more difficult measurement is to 
actually measure such a pulsar's proper motion. The direction 
of motion can confirm (or refute) the association with the SNR, 
while the magnitude of the motion gives an ind ependent esti- 
mate of the pulsar's age ( paensler & Frail 2000| ). 

PSR B 195 1+32 is a rapidly spinning (P = 39.5 ms) radio, 
X-ray, and 7-ray pulsar, l ocated on th e edge of the unusua l 
SNR CTB 80 (G69.0+2.7) dStrom 1987t iKulkarni et al. 1988b. 



While B 195 1+32 is spinning only slightly slower than the Crab 
pulsar, its low period derivative implies a much larger charac- 
teristic age, t c = P/2P = 107 kyr, and a much low er surface 
magnetic field, B«5x 10 11 G ([Fruchter et al. 1988b. 



and the large age for B 195 1+32 compared to most pulsar/SNR 
associations. 

While this evolutionary picture is aesthetically pleasing, it is 
imperative that it be verified by observations. Specifically, the 
model proposed for CTB 80 can be tested by measurement of 
the proper motion of PSR B 195 1+32. If the two sources are 
associated, then the pulsar should be moving away from the 
cent er of the SNR, as defined by the shell seen in infra-red and 
H I ( |Fesen, Shull, & Saken 1988| ; |Koo et al. 1990| ). The separa- 
tion between the pulsar and the shell's center, when combined 
with the pulsar's characteristic age, lets us predict a proper mo- 
tion for the pulsar of \x ~ 15 mas yr" 1 at a position angle (PA) 
w 250° (north through east), independent of the distance to the 
system. 

PSR B 195 1+32 shows significant timing noise, which pre- 
vents me asurement of its prop er motion through time-of-arrival 
analysis ( |Foster et al. 1994 ). Thus, only via interferometric 
measurements can this prediction be tested. In this Letter, we 
report on measurements of the motion of PSR B 195 1+32 over 
an 1 1 -year baseline. In §|2| we describe our observations, while 
in §|| we explain the method used to make our measurements 
and present our results. In §[| we discuss the implications of 
these results for the pulsar and SNR. 

2. OBSERVATIONS 

Observations of PSR B 195 1+32 were made with the Very 
Large Array (VLA) at epochs 1989.04, 1991.55, 1993.02 and 
2000.90. Each observation was eight hours in duration and used 
the VLA's A configuration. In the 1989, 1991, and 1993 obser- 
vations, data were recorded at two simultaneous frequencies, 
1385 & 1652 MHz, within the 20-cm band. Each frequency 
consisted of 15 channels spread across a 23.5 MHz bandwidth 
with a sampling time of 10 seconds, resulting in a useable field- 
of-view of radius ~ 10' centered on the pulsar. For the 2000- 
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The distance to the pulsar of 2.4 ± 0.2 kpc, as estimated from 
its dispersion measure, is consistent with the SNR's distance 
of 2 kpc as measured with H I absorption ( Strom & Stappers 
2000), suggesting that the two sources are physically associ- 
ated. However, this system appears quite different from more 
typical pulsar/SNR associations, in which a young (r c <; 20 kyr) 
pulsar sits near the center of an approximately circular shell. 
This difference can be understood if one invokes an evolution- 
ary scenario in which the pulsar was originally interior to the 
SNR, but as a result of its high space velocity has caught up 
with and begun to penetrate its SNR (Fesen, Shu ll, & Saken 
1988; |Hester & Kulkarni 1988| ; |Koo et al. 199$ . The pul- 
sar's relativistic wind then interacts with the SNR she ll, re- 
energizing and distorting it ( [Shull, Fesen, & Saken 1989 ). This 
sequence can explain both the strange appearance of CTB 80 
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epoch observations, we also incorporated data from the Pie 
Town antenna of the VLBA, which doubled the resolution of 
the array primarily in one dimension. The observations for that 
epoch were conducted at 1385 & 1516 MHz, and consisted of 
15 channels across a 12.5 MHz bandwidth with 5-sec sampling, 
again resulting in a 10-arcmin field-of-view. The flux density 
scale of our observations was determined using observations of 
3C 286, while the time-varying gains for each antenna were 
measured using regular observations of PKS B 1923+2 10 or 
TXS 2013+370. 

After standard editing and calibration, we produced images 
of the pulsar field for each epoch and frequency, using multi- 
frequency synthesis to mitigate bandwidth smearing, and dis- 
carding all baselines shorter than 10 km (corresponding to spa- 
tial scales larger than 4"). By removing these shorter baselines, 
we ensured that emission from SNR CTB 80 and from the com- 
pact wind-driven nebula surrounding the pulsar were not de- 
tected. The only emission seen in our images was the pulsa r 
itself, the "hot-spot" seen immediately adjacent (Strom 1987), 
and various other point sources spread throughout the field. 

We have identified seven such sources within 10' of the pul- 
sar, as listed in Table |l[ After deconvolving each image using 
the CLEAN algorithm, we applied a gaussian fit to the pulsar 
and to each of these seven sources to measure their position 
and extent at each epoch and frequency. 

3. MEASUREMENT OF THE PULSAR PROPER MOTION 

In order to accurately measure the pulsar proper motion, we 
adopted six of the field sources as reference sources, reserv- 
ing the nearest source to the pulsar (source 1 in Table [j]) as 
a check on our measurements. The reference sources are ap- 
proximately evenly distributed throughout the field. The qual- 
ity of each source as an astrometric reference was determined 
by measuring the vector separation between all possible pairs 
of sources, and then verifying that no source showed trends 
or significant changes in the magnitude and direction of these 
vectors. We also verified that in all cases the dimensions of 
the fitted gaussian matched those of the synthesized beam, con- 
firming that scatter broadening and bandwidth smearing (in the 
case of background sources) or scintillation (in the case of the 
pulsar) were not affecting our positional determinations. 

The ionosphere can potentially distort the measured posi- 
tions, producing a frequency-dependent angular displacement 
A9 = k6\ 2 , where 9 is the angular separation between a source 
and the field center, A is the observing wavelength, and k is a 
constant. By measuring 9 for each source and each frequency, 
we found that k had a value consistent with zero, demonstrating 
that any ionospheric effects were dominated by other uncer- 
tanties in our measurements. 

McGary et al. (200 1 ) used the VL A to measure pulsar proper 
motions an order of magnitude smaller than that expected here, 
and found a variety of other ways in which the positions of 
reference sources can be distorted, including relativistic effects 
due to the Earth's motion, errors introduced by the VLA corre- 
lator and additional empirical corrections which had no simple 
explanation. At the lower precision required here, we have ac- 
counted for all these effects simply by measuring the scatter in 
the position of each reference source between epochs. 

We began by determining the vector distances between all 
possible pairings of the six reference sources. The standard de- 
viations of the components of each vector in Right Ascension 



and Declination were computed for the four epochs; these were 
taken as the uncertainty for a particular pairing. This uncer- 
tainty was then decomposed between the two sources in that 
pair, the relative contributions of the two sources to the joint 
error being weighted inversely by each source's signal-to-noise 
ratio. We were thus able to derive an uncertainty in the position 
of each source for each pairing; these uncertainties were aver- 
aged over the five possible pairings to determine the best mea- 
surement of the true uncertainty in each coordinate for each 
source. This analysis was performed separately for observa- 
tions at 1385 MHz and at 1516/1652 MHz. A mean reference 
position for each epoch and frequency was then determined by 
averaging together the positions of the six reference sources, 
weighting the contribution from each source inversely by its 
distance from the field center. 

The separation between the pulsar's position and this mean 
reference position was then measured for each epoch and fre- 
quency, combining in quadrature the calculated uncertainties in 
the pulsar position with those determined for the reference posi- 
tion. The results of these measurements are plotted in Figure |], 
and show a clear motion of the pulsar to the south-west. The 
flux density of the pulsar was higher during the 1993 observa- 
tions than at other epochs, resulting in smaller uncertainties in 
the pulsar position for this measurement. Note that the position 
of the pulsar listed in Table jl] i s consistent with previous de- 
terminations ( foster et al. 1994 ), but with larger errors due to 
the systematic errors discussed above. Such uncertainties in the 
absolute astrometry do not affect our proper motion measure- 
ment, which has been determined based on astrometry relative 
to nearby background sources. 

As an independent test of our approach, we have similarly 
measured the proper motion for source 1 . As is shown in Fig- 
ure |lj no change in position is seen for this source, demonstrat- 
ing that the motion measured for the pulsar is real and that the 
uncertainties have been realistically assessed. 

By applying a weighted least squares fit to the pulsar's po- 
sition at each epoch, we find a proper motion at 1385 MHz of 
/z = 31 ±5 mas yr" 1 , at PA = 250° ±7°, and at 1516/1652 MHz 
of /j, = 26 ± 6 mas yr" 1 at PA = 240° ±10°. These two mea- 
surements are consistent with each other within their uncer- 
tainties. After taking into account apparent motion of the pul- 
sar at the level of 5 mas yr" 1 due to differential Galactic rota- 
tion, 7 the combination of our two measurements yields a mo- 
tion /i = 25 ±4 mas yr" 1 at PA = 252° ±7°. The corresponding 
projected velocity is V t = (240 ± 40)o?2 km s" 1 for a distance 
2d2 kpc. 

4. DISCUSSION 

The transverse velocity we have inferred for PSR B 195 1+32 
is in agreement with the value V, ~ 300 km s" 1 implied b y scin- 
tillation in its dynamic spectrum ( fouchter et al. 1988 ). Fur- 
thermore, as shown in Figure ^ the measured position angle of 
the proper motion agrees closely with the direction predicted 
if the adjacent "hot- spot" is interp reted as a bow-shock driven 
ahead of the pulsar (Strom 1987). It can also be seen in Fig- 
ure ^ that the pulsar is clearly traveling directly away from the 
center of its SNR. While this result is not unexpected given the 
already strong evidence for an association between the pulsar 
and SNR, it is surprising to realize that this is the first time such 
an effect has actually been observed. Such measurements for 
other young pulsars in SNRs have failed to demonstrate such 



To make this calculation, we have used the rotation curve of Fich, Blitz & Stark (1989) and have assumed a 10% uncertainty in the systemic velocity of the pulsar. 
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motion, and have either required specific m odels of SNR and 



pulsar evolution to main tain the association (Bailes et al. 1989 



Gaensler & Frail 200C ), or have caused the association to be 
abandoned (Thompson & Cordova 1994). 

While the direction in which the pulsar is moving is as pre- 
dicted, the magnitude of the pulsar's motion is not in agreement 
with the simplest expectations. By projecting the pulsar's posi- 
tion back by r c = 107 kyr to its inferred birth-site, it can be seen 
from Figure || that if the pulsar was born near the center of its 
SNR, then the pulsar's true age, t p , must be less than t c . 

We can make an initial estimate of t p by simply identifying 
the offset, 9, of the pulsar from the SNR's center. We estimate 
from the morphology of the infra-red shell shown in Figure g 
that O ~ 27', so that t p = 0//i ~ 64 kyr. However, a high veloc- 
ity progenitor and/or asymmetric SNR expansion can produce 
a significant offset betwe en the pulsar birth-sit e and the geo- 
metric center of the SNR (|R6zyczka et al. 1993|; D ohm-Palmer 
& Jones 1996; [Hnatyk & Petruk 1999j ). Given the complicated 
evolutionary history of CTB 80, it is perhaps overly optimistic 
to define a precise geometrical center and simply assume the 
pulsar was born there. 

If the birth-site is not at the SNR's geometric center, there is 
no reas on to expect the pul sar to be traveling away from this po- 
sition ( |Gvaramadze 2000 ). However, the farther the birth-site 
from the center of the SNR, the less likely the pulsar should be 
moving in this direction by chance. We can therefore quantify 
the uncertainty in the pulsar's birthplace as follows. We assume 
that for a given angular separation <j> from the SNR's geometric 
center, it is equally likely that the pulsar was born at any posi- 
tion on a circle concentric with the SNR and of radius <fi. Thus, 
for a particular value of <f>, we can compute the fraction of pos- 
sible birth-sites which are consistent with the pulsar's direction 
of proper motion, PA = 252° ±7°. By considering the full range 
of values <(j>< 8, we can build up a probability distribution 
of possible pulsar ages, from the integral of which we can de- 
termine the lcr confidence limits on <j>. Using this approach, we 
find that the pulsar was born within ±8' of the SNR's center, 
and that the age implied by our proper motion measurement is 
t p = 64 ± 18 kyr. These estimates are independent of the dis- 
tance to the system and on the site of the supernova explosion; 
they depend only on the assumptions that an approximate ge- 
ometric center can be defined for the SNR, and that the pulsar 
is moving away from its birthplace. The age we have derived 
is consisent with the estimate of 77 kyr, derived from the ex- 
pansion v elocity of the H I shell coincident with the SNR (K oo 
et al. 1990). 

It is not surprising that this system's true age is less than the 



pulsar's characteristic age. For a pulsar of constant moment 
of inertia a nd magnetic dipole momen t, it is straightforward to 
show that ( [Manchester & Taylor 1977 ): 



t B = 



2t c 



n-\ 



(1) 



where n is the pulsar's "braking index" and Pq is its initial spin- 
period. For n = 3 and Pq <C P, Equation (|l|) reduces to t p = t c , 
but in general neither such condition will be satisfied. While n 
has not been measured for PSR B 1 95 1 +32, we find 1 .5 <; n <; 3 
for the five pulsars which have had their braking indices mea- 
sured (Zhang et al. 2001 and references therein). For n = 3, 



the range in t p determined above implies Pq = 25 ± 5 ms for 
PSR B 195 1+32, while for n = 1 .5 we find P = 29 ± 3 ms. Com- 
bining these estimates, we find the initial period most likely 
falls in the range Pq = 27 ± 6 ms. Again, this estimate does not 
depend on the assumed distance to the source. 

As listed in Table g, PSR B 195 1+32 joins a rapidly grow- 
ing list of pulsars whose initial periods have been determined 
from their associations with either SNRs or historical super- 
novae. While these measurements span a reasonable range of 
initial periods, the data provide good evidence that radio pul- 
sars are born as rapid rotators, with no evidence for a popula- 
t ion of longer initial periods a s has been sometimes propose d 
( [Vivekanand & Narayan 1981 ; Chevalier & Emmering 1986). 
On the other hand, it is becoming abundantly clear that r c can 
be an unreliable estimator of a young pulsar's age, with corre- 
sponding implications for associations with SNRs, pulsar ve- 
locities, and models for neutron star cooling. With many new 
pulsar/SNR associations now being identified at both radio and 
X-ray wavelengths, a statistically useful sample of pulsar initial 
spin periods should soon emerge. 
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Table 1 

Positions of sources at epoch 2000.90. Values in parentheses indicate the uncertainty in the last digit. 



Source 


Position (J2000) 


Distance from 




Right Ascension 


Declination 


pulsar (arcmin) 


PSRB1951+32 


19 h 52 m 58?2(l) 


+32°52'41"(1) 




1 


19 h 53 m 17?80(l) 


+32°52'09."6(1) 


4.1 


2 


19 h 53 m 16?36(l) 


+32°48'46"4(1) 


5.5 


3 


19 h 53 m 15?63(l) 


+32°59'38"1(1) 


7.9 


4 


19 h 53 m 25?41(l) 


+32°58'12"5(1) 


8.0 


5 


19 h 53 m 13: s 41(l) 


+33°01'21."8(1) 


9.3 


6 


19 h 52 m 15?79(l) 


+32°49'35"7(1) 


9.4 


7 


19 h 52 m 15?64(l) 


+32°49'39"2(1) 


9.4 
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Table 2 

Initial period estimates for young pulsars. 





rYSSOClalcU olM\ 


lillLlal r CI1UU 


JVClClCllOC 




or Nebula 


(ms) 




J0537-6910 


N 157B 


< 14 


Marshall et al. (1998) 


B0531+21 


Crab Nebula 


19 


Manchester & Taylor (1977) 


B1951+32 


CTB 80 


27 ±6 


This paper 


B0540-69 


0540-69.3 


30±8 


Reynolds (1985) 


J0205+6449 


3C58 


60 


Murray et al. (2002) 


J1811-1925 


G 11.2-0.3 


62 


Toriietal. (1999); Kaspi et al. (2001) 


Jl 124-5916 


G292.0+1.8 


90 


Camiloetal. (2002) 
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FIG. 1 . — Proper motion measurements of PSR B 195 1+32. The data points show the offset of the pulsar from the mean reference position at four epochs, and 
the lines show the weighted best fit to these data. The solid circles and solid line correspond to 1385-MHz data, while the open circles and dashed line represent 
1516/1652-MHz data. From left to right, the data points correspond to observations at epochs 1989.04, 1991.55, 1993.02 and 2000.90. The inset shows the position 
of source 1 at each epoch, plotted on the same scale and determined via the same process. 




RA (J2000) RA (J3000) 

FIG. 2.— SNRCTB 80 and PSR B 195 1+32. The left panel presents an image of the IRAS 60/100 fim ratio in the region, showing the shell of CTB 80. The current 
position of PSR B1951+32 is marked by the "+" symbol, the direction along which it has traveled is indicated by the solid line, and its projected birthplace for an 
age of 107 kyr is marked by the lcr ellipse. In the right panel, we show a 1.5-GHz VLA image of the immediate vicinity of the pulsar, made from all baselines of 
our 1993-epoch data. PSR B1951+32 is the point-source to the south-west of center. The arrow indicates the pulsar's measured direction of motion, the length of 
the arrow corresponding to the distance traveled by the pulsar over 100 years. 



